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Abstract 

Graviton production due to collapsing extra dimensions is studied. The momenta lying in the 
extra dimensions are taken into account. A D-dimensional background is matched to an effectively 
four-dimensional standard radiation dominated universe. Using observational constraints on the 
present gravitational wave spectrum, a bound on the maximal temperature at the beginning of the 
radiation era is derived. This expression depends on the number of extra dimensions, as well as 
on the D-dimensional Planck mass. Furthermore, it is found that the extra dimensions have to be 
large. 



1 Introduction 

Our observable universe is four-dimensional. However fundamental theories, like for example string 
theories, require for consistency more than four space-time dimensions. In the standard Kaluza- 
Klein approach EJ the extra dimensions are compact and curled up to very small size, roughly of the 
order of the Planck scale, £p ~ 10 -33 cm, implying that space-time is effectively four-dimensional. 
However recently, with the discovery of brane solutions to string/M-theory, the idea of large extra 
dimensions has attracted much attention |||. In this last approach, it is proposed that ordinary 
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matter is confined to a three-dimensional sub-manifold, a three-brane, which is embedded in a higher 
dimensional space-time; the graviton on the other hand is allowed to propagate freely through the 
whole space-time. Neglecting the brane tension, i.e. the energy density per unit three-volume of the 
brane, and considering compact dimensions, one re-introduces the Kaluza-Klein picture. However, 
in this case the extra dimensions do not have to be small. Newtonian gravity has been tested and 
found to hold down to scales of order of 1 mm @. Below this scale, gravity could in principle 
be higher than four-dimensional. In higher dimensional gravity, the four-dimensional Planck scale 
is no longer a fundamental scale; the higher dimensional Planck scale, Md, becomes instead the 
fundamental scale. This allows to explain the huge difference between the electroweak and the 
four-dimensional Planck scale, known as the hierarchy problem. Assuming for simplicity that the 
n extra dimensions form an n-torus which has the same radius R in each direction, and using 
Gauss' law, the D-dimensional and the four-dimensional Planck masses, M4 and Md respectively, 
are related by (§ 

Ml = R n M™ +2 , (1.1) 

where R n is the volume of the n extra dimensions. Thus, taking Md to be of the order of the 
electroweak scale, Md ~ ITeV, the huge difference between M4 and the electroweak scale can be 
explained as due to the large size of the extra dimensions. 

Here we are going to study the effects of time- varying extra dimensions on the graviton produc- 
tion. We assume that the background space-time can be written as a direct product of an internal 
space, i.e. the extra dimensions, and a four-dimensional external space-time. Furthermore, during 
an initial phase, the internal space is contracting and the four-dimensional external space-time 
is expanding. At some time this is matched to a radiation dominated flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe with the internal dimensions frozen at constant size. Due to 
the changing background metric, there will be particle production. Here we will be in particular 
concerned with the production of gravitons. Using observational bounds on graviton spectra it is 
possible to derive a relation between the D-dimensional Planck mass and the temperature at the 
time of transition. Usually it is assumed that there are no massive modes excited, i.e. momenta 
lying in the internal space are not taken into account. Here however, we are going to assume that 
the internal momenta are excited and we are interested in their effect on the spectral energy density 
in four dimensions. In order to determine this effect, the internal momenta are integrated out. This 
will lead to a final expression for the spectral energy density, which will be just depending on the 
four-dimensional momenta. 

The formalism to deal with metric perturbations has already been developed in Ref. || . Particle 
production in higher dimensional space-times has been also the subject of Refs. [f?J and ||. 

2 Graviton production 

Consider a D{= d + n + l)-dimensional space-time with line element, 

ds 2 = a 2 ( V ) [drf - JydsW] - b 2 {r,)8 AB dy A dy B , (2.2) 

where d = 3, and the indices i, j = 1, .., 3 and A, B = 4, .., 3 + n. Thus, a and b denote the scale 
factors for the internal and external spaces respectively, 77 stands for the conformal time, and 5ij 
is the Kronecker symbol. Assume that the extra dimensions are compact and, for simplicity, that 
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they all have the same size, i.e. < y < 2-kt. Then the length-scale, R, characterizing the physical 
size of the extra dimensions is R = lixrb. 

The higher dimensional phase is matched to a radiation dominated universe. It is assumed that 
the extra dimensions are dynamical during the first stage. At the onset of the radiation dominated 
era they become static. 

Hence, the evolution of the background will be parametrized as follows, 

b(v)=h[-— ) , for 7/ < -r]i 

b{rj) = b\ , for r\ > —r)\ . (2-3) 

In the following we set a\ = b\ = 1. In Figure 1 we present the evolution of the background for 
n = 8 extra dimensions. 
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Figure 1: The evolution of the background during a twelve-dimensional and an effectively four- 
dimensional stage. Here the transition takes place at rj = —771 = —1., and the number of extra 
dimensions is given by n = 8. The solid line shows the scale factor a(rj) of the four-dimensional 
space-time and the dashed line shows the scale factor b(rj) of the eight-dimensional extra space. 

Within the context of the Kaluza-Klein theory, there are expected oscillations of the 6-field 
about the minimum. This is not considered here, since this issue represents a question on its own. 
As it can be seen from Eq. ( |2.3[) , we consider a power-law behavior for the 6-field. Parametric 
resonance in the context of a higher Z)-dimensional Kaluza-Klein theory has been only studied, to 
our knowledge, for Kaluza-Klein modes ||. 

The tensor perturbations of the D-dimensional metric G ^ take the form || 

/ \ 

6G AB = -a 2 7ii(^,y) 

V -b 2 lAB (x,y) J 

We consider that the unperturbed D-dimensional metric is a direct product, thus there are no cross- 
terms in the above expression for the tensor perturbations. Only 7^ describes gravitational waves 
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in the four-dimensional space-time. Although ^ab describes gravitational waves in the n extra 
dimensions, it is interpreted as a scalar (matter) perturbation in the four-dimensional space-time 
of the observer. Both types of gauge-invariant tensor modes satisfy the D-dimensional Klein-Gordon 
equation for a massless scalar field. Here, we are only interested in the modes which behave as 
gravitational waves in the four-dimensional space-time. Writing 7^ = (j)(t,x,y)eij, where is the 
polarization tensor, the amplitude of the gravitational waves <f> satisfies the same equation as a 
minimally coupled scalar field, i.e. || [ 10 1 



□ = , 

where □ denotes the d'Alembert operator in D-dimensions. In the above background, in Fourier 
space, this yields to the mode equation for the canonical field <I> = ab n / 2 <j), 



k 2 + 
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*Z = , 



(2.4) 



where the three-vector k lies in the external space, the n-vector q is in the internal (compact) 
space, I stands for (k, q), and 

j3 = a-X 

(ab n / 2 )" V 2 



N 



a (a - 1) + naX + n(n - 2)A 2 /4 + nX(X - l)/2 . 



ab n / 2 

Furthermore, the modulus of the external and internal momentum vectors is denoted by k = \/|k| 



and q = \/|q[, respectively. Using the compactness of the extra dimensions the components q A can 
be written as q A = Na/^, where Ma are integers. This leads for the modulus of q to 



M 
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(2.5) 



with J\f given by J\f = yYl-^A- ^ or @ = — ^ exact solutions can be found, which are discussed in 
Ref. ||. However, for general values of (5 there is no known exact solution. Therefore, we follow 
Ref. 1 11] to find an approximate solution. In the following, it is assumed that a < and A > 0, 
which implies (3 < 0. This describes a universe with a collapsing internal (extra) dimensions and 
expanding external dimensions, before the transition to an effectively four-dimensional space-time. 



2.1 Mode evolution in the presence of dynamical internal dimensions 

During the higher dimensional state the background is given by the vacuum Kasner metric. This 
constrains the exponents A and a, and thus /3, and allows to write the metric as a function of the 
number of extra dimensions n. 

Namely, the Kasner conditions da^ + na\ = 1 and da\ + naf = 1 imply for d = 3 and 
that 
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n 



3 + n 



-n)< 



'3 2 + n 
n (3 + n) 2 



±J 3 - 2 + n 



3 + n y n (3 + n) 5 



(2.6) 
(2.7) 
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These are related to the exponents a and A by 



1 - «E 
1 — «E 



A = T ^_. (2.9) 



Hence, calculating the parameter N in Eq. (2.4) yields N = —(1/4), independently of the number of 
extra dimensions. This can be understood by realizing that ab n l 2 = (a 3 6 n /a) 1//2 = (Volume/a) 1//2 , 
where the volume is constrained by the first Kasner condition. 

The parameter (3, which describes the difference in the expansion rates of the internal and 
external space, takes values —1 < (3 < —1/(1 + v3), where the lower limit corresponds to n = 1 
and the upper one to the limit of a large number of extra dimensions n. 

The approximate solution of the mode equation, Eq. (|2.4j ), is found by solving it in two different 
regimes. Following the same approximation as in Ref. [O], two cases emerge. 



Case (i): Assume that in Eq. ( |2.4[) the square of the three-momentum k always dominates over 
the term involving the internal momentum q. For very early times, r] — > — oo, or equivalently 
very large k, hence very short wavelengths, k 2 also dominates over the ?7 _2 -term. And thus, the 
solutions are (in-coming) plane waves. However, for times closer to the transition at rj = —rji, or 
for smaller k, hence longer wavelengths, k 2 < 1/rj 2 . Hence, this implies for the modulus of the 
internal momentum q, that q < (— 77/771 ) _ ^/c, and thus together with rj\ = l/fci, the condition 
reads: q/k < {k/k\)". Note that this case also includes the case q = 0. The in-coming state is 
effectively four-dimensional. Then Eq. ( |2.4D reduces to a Bessel equation which is adequately, i.e. 
correctly normalized to an in-coming vacuum state, solved by 



*z = Vt^£ 2) (IN) for r i<-m, (2.io) 

with the Hankel function of the second kind and /i 2 = (1/4) + N. Together with A^ = —1/4 
the index of the Hankel function is identically vanishing, fj, = 0. This solution is matched on super- 
horizon scales, i.e. krji <C 1, to a radiation dominated universe at r] = —rji. Thus, the canonical 
field $z and its first derivative are given respectively by 

2 

®i\ v =-m ~ -i-Vmhi(krji) 

*'i\n=-m ~ --Vihl + HkVi)] ■ (2.H) 

7T 



Case (ii): Assume that the q-term becomes dominant before the perturbation in D dimensions 
becomes super-horizon, hence k 2 + q 2 (—rj/r]i) 2 ^ > l/r? 2 . Thus, the mode equation can be approxi- 
mated by 



®i + 



k'+qH-^ 

m 



23 



$i = o, 



(2.12) 
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which is approximately solved by (cf ]Tl[) 



exp 



irj\ k 2 + K 2 q 2 
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(2.13) 



where k = 1/(1 + /3) , j3 / —1. The case (3 = — 1 can be solved exactly and it is discussed in Ref. Q. 

At conformal time r? = i] q the internal momentum becomes dominant over the external mo- 
mentum and finally the rj~ 2 -tevm dominates. For 77 > rj q the mode equation can be approximated 
by 



m 



2/J 



which is solved by fill 



N 

7^2 



(2.14) 
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\qrj\K 



gr/ «; 



(2.15) 



Matching Eq. ( [2.15 ) to Eq. ( 2.13j ) leads to c q L> = and Cq = ij^fq. Hence, for \qi]i\ <C 1, i.e. 

(2) 

super-horizon perturbations, and using 

Hq (*) i(2/vr)ln z, one obtains that and its first 

derivative at r\ = — rji are given respectively by 

<&l \r)=-r)i ~ 



$z l«7=-m ~ 



2 i 
■i—rji In(Kqrji) 

7T 

i -i 2 -i 

— 77, 2 ln(KCTi) — i— (/3 + 1)771 2 . 

7T 7T 



(2.16) 



2.2 Mode evolution with static internal dimensions 

The higher dimensional state evolves instantaneously into an effectively four-dimensional radiation 
dominated universe at 77 = —771. The extra dimensions stay at a constant size after the transition, 
with scale factors given by Eq. (|2,3] ) for 7/ > —771. Then the mode equation, Eq. (|2.4|), takes the 
form 



®l + 



$,. = 0. 



(2.17) 



This can be transformed into the equation for parabolic cylinder functions by introducing z = 
(217/771 J 1 / 2 (77 + 2771) and a = -(rn/2q)k 2 |]. This yields 



d 2 ®i 
dz 2 



+ 



a = 0. 



Its complex standard solution is given by |1J 

$ l (z)=B[c-E(a,z) + c + E*(a,z)} , 



(2.18) 



(2.19) 
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where for later convenience the constants are chosen to be the Bogoliubov coefficients c± and B 
is a normalization factor. With the Wronskian condition on the mode functions (see for example 
Ref. |5J), i.e. $/— = i, and the normalization of the Bogoliubov coefficients |c + | 2 — |c_| 2 = 1, 
together with the expression for the Wronskian of E(a,z), E*(a,z) JT^], the normalization factor 
B is found to be B = [r ]l /(2q)} 1 / 4 /y/2. 
Hence the mode solution is given by 



*(*) 



V2 \2q 



[c-E(a,z) + c + E*(a,z)] 



(2.20) 



In order to make progress the solutions will be considered in two regimes of approximation which 
will comprise cases (iii) and (iv). 



Case (iii): Assume that \a\ > z 2 /4 44> k 2 > [(rj + 2rji)/rji] 2 q 2 . Using the limiting behaviour 
given in Ref. Jl2| ([19.22.4]), E can be approximated in this limit as 



E{a,rj) ~ y/2 



2 3 / 4 



ra + Ua) 



r(4 + Ua) 



exp 



1 / q\ 2 I r] + 2r/i 



4 \k 



Jk(ri+2rn) 



(2.21) 



and a similar expression for the complex conjugate E*. Using 



r(| + Ua) 



r(f + \ia) 



the mode functions in this limit evolve according to 
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(2.22) 



Case (iv): Assume that |a| < z 2 /4 44> fc 2 < [(77 + 2ryi)/ryi] 2 /g 2 . With the approximation of 



E(a, z) |12[ in this limit the evolution of the modes is given by 



2q(r 1 + 2ri l ) 



C- exp 



z J_( r? + 2??1 ) 2 
2171 



+ c_|_ exp 



-i-^-(r / + 2r/ 1 ) 2 
2771 



(2.23) 



2.3 Graviton spectra 

Gravitons are produced due to the changing background space-time. Their spectral energy density 
might allow conclusions about the early history of the universe. The Bogoliubov coefficients trans- 
form the in-state into the out-state and determine the number of created particles. As it turns 
out there are three possible combinations of the four mode solutions discussed previously. It is 
found that the evolution of modes from the regime of validity of cases (i) and (iii), (i) and (iv), 
as well as (ii) and (iv) is possible. For these three regimes, the Bogoliubov coefficients have to be 
calculated. Then, the spectral energy density is found by simply adding up the contributions from 
the three parts. In order to find the total spectral energy density in the effectively four-dimensional 
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space-time after the transition, the momenta lying in the internal (extra) space are integrated out. 
In this way the contribution of energy due to the modes in the extra space can be estimated. 

The Bogoliubov coefficients are calculated by matching the mode functions and their first deriva- 
tives at the transition time 77 = — 771, on super-horizon scales. The resulting expressions together 
with the appropriate range of validity are given below. 



Case (i) matched to (iii): 



This holds for | < a 1 . 



\/2tt \/kr)i 



Case (i) matched to (iv): 



1 Ml) 2 
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(I) +2ik m ) ln(k Vl ) 



(2.24) 



■ QV1 

c_e l ~ r __ [l + (l + i q7]1 ) ln(feT?i)] , 



(2.25) 



which holds for a 1 < f < (^j ^ . 



Case (ii) matched to (iv): 



c_e 2 11 ~ — 



V2_ 



[(3+l + (l+iqrn)ln(Kq m )} 



(2.26) 



This holds for § > (fy' 



Case (ii) matched to (iii): 

Clearly, this possibility is not allowed within our set-up, where the extra (internal) dimensions are 
collapsing while the external four dimensions are expanding. This implies (3 < 0, and therefore 
case (ii) requires q/k > 1, while case (iii) holds for the opposite limit, namely qjh < 1. Therefore, 
our scenario does not allow the matching of cases (ii) and (iii). 



In order to compare with observations it is useful to calculate the spectral energy density. The 
total energy density of created particles is given by 



(2< 



n+3 



a? + b 2 



dV r 
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(2.27) 



and dVphys = dV com / (a s b n ), where it is assumed that the comoving volume consists of two spheres, 
i.e. 
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which for d = 3 implies 
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(2.28) 



Hence the spectral energy density p(k) = dp/dlogk is given by 

4 / t \ n 



p(k) 



R n 



-16 



i+- 

7T + 2 



feV fk 



(2^+3 ^rm V«y U 



dry 



n-l 



1+1 1. 1 FJ 



(2.29) 



where Y = q/k. 

As explained above there are three contributions to the spectral energy density. The notation 
P(i)—(iii) stands for the contribution resulting from the Bogoliubov coefficients calculated from the 
matching of cases (i) and (iii), and in a similar way for the other cases. For pu^—uu), the integral 
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has to be calculated. For pu)-a v ), the corresponding one. In order to calculate P(u)-(i V ), we 
introduce an upper cut-off in the integral, i.e. 
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At the transition time r\ = —r\\ the metric is continuous but its first derivative is not. This leads 
to the so-called sudden transition approximation. It basically means that for the modes with 
periods much greater than the duration of the transition phase, the transition can be considered 
as instantaneous. It is known that this type of approximation leads to an ultra-violet divergency 
0jn|. This justifies the introduction of an upper cut-off q mSLX . Neglecting sub-leading terms for 
the total spectral energy density, p(k) = po^-im) + P(i)-(i v ) + P(u)-(iv)i the following expression is 
found 
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(2.30) 



This can be expressed in terms of the fractional energy density (dimensionless parameter), Q,(k) = 
p(k)/p c , where the critical energy density p c is given in terms of the D-dimensional Planck mass, 
Md, as p c = [3/(8vr)]# n M™ +2 # 2 . We thus find that the fractional energy density of produced 
gravitons OgwC 7 ?) a t a given time rj, with respect to the fraction of critical energy density in 
radiation, U y (r]), defined as ^ 7 (f?) = (Hi / H) 2 (a\ / a) 4 , is 



itawW ~ Aa-«(^) 2+ "n 7 (,,(*) 3+ "(m* 
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+ _4 2 „(__) (ln- 

+ ■ 4 -(^) 2+ " f! >»(i) S (^)"( taK 'f £ ) 2 - ( ^ 31) 

The constants Ai contain all the numerical factors. Furthermore, it was assumed that the maximal 
frequency k\ in the four-dimensional space-time is of the same order as the curvature scale at the 
time of transition, i.e. k\ ~ H\. 

This expression can be compared with the one for no excited internal momenta. It can be easily 
recovered from the matching of case (i) to (iii), without carrying out the integration. It is found 
that 

«cw , =A[§-) 2 H, (,) (£) 3 (in |) 2 . (2.32) 



As it can be seen from the expression for the fractional energy density, Eq. ( 2.31| ), the effect of 



the momenta in the extra dimensions is two-fold. On the one hand, there are two contributions 
which are growing with the scale factor, and therefore might soon dominate the energy density of 
the universe and overdose it. On the other hand, the cut-off g max has to be chosen in such a way 
that its contribution does not dominate, i.e. Qmax ^ (/^ -Q^i ■ Thus, the cut-off of the momenta 
in the extra space has to be less than the one in the four-dimensional space-time, q ma , x < 

There are several constraints on the gravitational wave spectrum. The strongest ones come 
from big bang nucleosynthesis and the prevention of recollapse of the universe. Both of these hold 
for the whole range of frequencies. Big bang nucleosynthesis successfully predicts the production 
of light elements. However, everything depends crucially on the expansion rate at the time of 
nucleosynthesis. This in turn is determined by the total energy density at that time. Therefore, 
any extra contributions at that time are constrained. The total energy density of gravitational 



waves at the present time has to satisfy [[15| [|lq 1 

/ d(log/) h 2 n G w(/) < 0.227 Q 7 , (2.33) 

Jf=Q 

where 0.5 < ho < 0.85 parametrizes the experimental uncertainty for the present value of the 
Hubble parameter Hq. One might argue that, as a first approximation, /iq ^gw(/) ^ 0.227 f2 7 , 
which leads to 

/ tt \ 2+n 

A 2 a < 0.227 . (2.34) 



Furthermore, ao is given by 



H,\ 112 _ ^„/H^ 



where the subscript denotes the present epoch. Assuming that the standard behaviour of radiation 
holds up to the transition scale, the Hubble parameter H and the temperature T are related by 
(see for example p7[|) 

H = lMgl^, (2.35) 
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Figure 2: The maximal temperature at the beginning of the radiation dominated FLRW stage as 
a function of the D-dimensional Planck mass. The solid line shows n = 2 extra dimensions, the 
long dashed one is for n = 6, the dashed-dotted line for n = 7 and the long/short-dashed one for 
n = 10. 



where g*(T) counts the total number of effectively massless degrees of freedom. For temperatures 
T >300 GeV, one has = 106.75. Using Eq. ( 2.35|) for the value of the Hubble parameter at the 
time of transition, H±, the bound given by Eq. fl2.34| ) can be written in terms of the temperature 
at the transition, i.e. at the beginning of the radiation dominated FLRW stage. The maximal 
temperature T\ is given when the bound is satisfied. It is found to be 
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(2.36) 



where / = 0.227 is our bound. This is shown in Figure 2. It should be noted that the range of the 
D-dimensional Planck mass is constrained by the upper size of the extra dimensions, i? max ~ 1mm. 
This gives the constraint 



M D > (2 x lO^ 32 /™)- 16 



V(i+S) 



TeV 



(2.37) 



Therefore, there is a lower bound on T\. The minimal value of M-q and the minimal value of T\ 
as a function of the number n of extra dimensions are shown in Figure 3. Assuming that T\ is 
the same as the reheating temperature Trhj Eqs. ( p. 36 ) and ( |2.37| ) imply a lower bound on Trh- 
More precisely, for one extra dimension (n = 1), the lower bound for the L>-dimensional Planck 
mass is Md > 10 5 TeV, leading to a lower bound for the reheating temperature Trh > 3 x 10 9 GeV. 
In a similar way, for two extra dimensions (n = 2), one obtains M D > ITeV and T RH > 10 7 GeV. 
Finally, for n = 3 one has M D > 10~ 3 TeV and T RH > 5 x 10 6 GeV. 

Big-bang nucleosynthesis must proceed in the standard way, implying that the reheating tem- 
perature must be higher than ~ lOMeV. This requirement is satisfied in our model, as the above 
analysis shows. Standard inflationary models, on the other hand, lead to a higher reheating tem- 
perature than the lower bounds we found above. However, one must keep in mind that above we 
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Figure 3: In the left panel the allowed minimal value of the D-dimensional Planck mass Md is 
shown as a function of the number of extra dimensions n. In the right panel the minimal value of 
the temperature at the beginning of the radiation dominated stage, T\, is shown as a function of 
n. 



only gave the lower bounds, while the actual reheating temperature can be much higher indeed. In 
addition, a higher value can be easily achieved increasing the value of the D-dimensional Planck 



mass, provided it always satisfies the constraint given in Eq. ( 2.37 ). Finally, reheating from a 
higher number of collapsing extra dimensions favours larger values of Mq, and thus smaller extra 
dimensions. 

The other strong constraint comes from the requirement that the additional energy density 
should not overdose the universe. Thus, the condition £Igw(lu) < 1 has to hold for all frequencies. 
However, if the big bang nucleosynthesis constraint is satisfied, then indeed Qqw < 1- 

There are additional constraints coming from large scales as probed by the cosmic microwave 
background explorer (COBE). Strong perturbations can change the gravitational potential experi- 
enced by the microwave background photons. This causes a red-shift of these photons and hence 



a fluctuation in temperature of the cosmic microwave background [15]. There is, in addition, a 
bound due to pulsars whose signal is delayed, if a gravitational wave passes through between the 
earth and the pulsar [^] . However, both these constraints apply only for a range of frequencies and 
they are in general less tight than the ones mentioned earlier. Furthermore, since the spectrum is 
increasing with increasing frequency, the most stringent bounds will come from large frequencies, 
i.e. small wavelengths. 

In conclusion, we find that the excited internal modes are compatible with observations, pro- 
vided the temperature at the transition time to an effectively four-dimensional universe, is less than 
its maximum, given in Eq. ( |2.36| ). Thus, recalling that Hi ~ T^/M^, H\ ~ k±, g m ax < ki and 
q ~ -R" 1 , this implies that the size R of the extra dimensions has to be R > M^/T^. As we have 
shown, the lower bound of the reheating temperature depends on the number n of extra dimensions 
and the value of the D-dimensional Planck mass Md . With values of Md , such that the constraint 
Eq. ( |2.37| ) is satisfied, one can easily check that our analysis leads to the requirement that the size 
of the internal space must be much larger than the four-dimensional Planck length 

R > Ipi , (2.38) 
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where £p\ is the four-dimensional Planck length (£p\ ~ 10 GeV ). 

Our model could very well be realized within the context of the pre-big-bang scenario [14], a 
particular model of inflation inspired by the duality properties of string theory. In the pre-big- 
bang scenario, perturbations of Kalb-Ramond axions can provide a quasi-scale-invariant (Harrison- 
Zel'dovich) spectrum of large angular scale cosmic microwave background (CMB) anisotropies flT8| . 
An extension of this mechanism (for massless axions) to the region of the acoustic peaks, showed 
that a consistency with the current CMB data requires that the internal dimensions contract at 
a rate faster than the rate at which the external dimensions expand [19|. Therefore, combining 
our model with the one where the universal axion of string theory triggered the CMB anisotropies, 
we get further constraints on the extra dimensions of string theory. However, this will be left for 
future work. 



3 Conclusions 

We have discussed graviton production in the context of a higher dimensional model, which from a 
multi-dimensional phase where the extra dimensions are contracting and the external dimensions 
are expanding, enters into an effectively four-dimensional universe with static extra dimensions. 
The momenta in the extra space were taken into account. To find an estimate of their contribution 
to the energy density in the four-dimensional space-time, the momenta in the internal space were 
integrated out. 

The contributions of the internal momenta tend to dominate the gravitational wave spectrum, 
up to the point of forcing the universe to recollapse at an early stage. However, one can assume 
that there is an upper cut-off for the internal momenta, given by the maximal frequency of the four- 
dimensional (external) space-time. The gravitational wave spectrum is constrained by observations. 
With the assumption that the standard relation between the Hubble parameter and the temperature 
holds up to the transition time, we derived an upper bound on the temperature at the beginning 
of the radiation dominated era. Using that Newtonian gravity has been tested sucessfully down to 
scales of the order of 1 mm, one obtains a lower bound on the temperature at the transition scale, 
which depends on the number of extra dimensions. This lower bound on the reheating temperature 
is, in general, lower than the reheating temperature of standard inflationary models. However, it 
can be raised if one assumes smaller extra dimensions, which lead to a higher /^-dimensional Planck 
mass. Our final conclusion is that the size R of the internal space must be much larger than Planck 
length £p\, namely R^> £p\. 
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